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Csj ■ Abstract 

We consider the effect of discretization errors on the microscopic spectrum of the Wilson Dirac 
operator using both chiral Perturbation Theory and chiral Random Matrix Theory. A graded 



^ ' chiral Lagrangian is used to evaluate the microscopic spectral density of the Hermitian Wilson 

Dirac operator as well as the distribution of the chirality over the real eigenvalues of the Wilson 
Dirac operator. It is shown that a chiral Random Matrix Theory for the Wilson Dirac operator 
reproduces the leading zero-momentum terms of Wilson chiral Perturbation Theory. All results 
are obtained for fixed index of the Wilson Dirac operator. The low-energy constants of Wilson 
chiral Perturbation theory are shown to be constrained by the Hermiticity properties of the Wilson 
Dirac operator. 
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I. INTRODUCTION 



In the continuum, low-lying spectra of the Dirac operator for theories with spontaneous 
chiral symmetry breaking have two equivalent descriptions. One is in terms of chiral Random 
Matrix Theory jl-4|, and the other one is in terms of a chiral Lagrangian 5-7|. It is by 
now well-established how the two formulations are in one-to-one correspondence {g, 8, 9|, 
and that this is valid for all spectral correlation functions of the Dirac operator. Even single 
eigenvalue distributions can be derived in both formalisms, and have been shown to be 



equivalent 



10| . The equivalence is valid to leading order in a chiral counting scheme known 



as the e-regime or, in Random Matrix Theory terminology, the microscopic domain. 

Apart from its conceptual value, the theory of low-lying Dirac operator spectra has been 
of quite practical use in lattice gauge theory. In fact, it serves multiple purposes, all relying 
on a QCD partition function that is formulated from the outset at finite four-volume V: 
(1) it can be used to establish spontaneous chiral symmetry breaking in a clean way, (2) it 
provides precise non-perturbative analytical predictions that can be used to test the chiral 
limit, and (3) it allows for a determination of low-energy constants by means of finite- volume 
scaling. So far these uses have been limited by the fact that violations due to finite lattice 
spacings, a, have been ignored. Such lattice artifacts evidently depend on the particular 
lattice discretization chosen. Following Symanzik's program, corrections due to finite lattice 
spacings can, when they are sufficiently small, be analyzed in a continuum field theory 
language through the introduction of higher- dimensional operators in the Lagrangian. The 
corrections to the chiral Lagrangian that arise up to and including order a^-effects for Wilson 



fermions have been analyzed in a series of papers 



llMlq. For comprehensive reviews of 



effective field theory methods at finite lattice spacings, see, e.g., ref. [ij 
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It is then an obvious problem to investigate the effect of lattice-induced scaling violations 



on the spectral properties of the Wilson Dirac operator [19|. In a recent Letter 20|, 
three of the present authors have taken up this issue and shown how the microscopic scaling 
regime can be phrased both in terms of the (graded) chiral Lagrangian and a modified chiral 
Random Matrix Theory that incorporates finite lattice spacing effects of Wilson fermions. 

Although the Wilson Dirac operator is not Hermitian, the operator = '-/^{Dw + fn) 
is Hermitian, and much more convenient to work with in lattice QCD simulations. In this 
paper we analyze the microscopic spectrum of this operator. Contrary to the lattice QCD 
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Dirac operator at a = 0, its eigenvalues A|(m) are not paired and are nontrivial functions 
of the quark mass. 

There is a deep relation between the topology of gauge field configurations and the spec- 
trum of the Dirac operator. Not only is the number zero eigenvalues equal to the difference 
of the number of right handed and left handed zero modes, because of level repulsion the 
Dirac spectrum near zero is affected in a universal way by the topological charge. As can be 
shown from spectral flow arguments (see, e.g. 2lN27j]). the eigenvalues corresponding to the 
chiral modes at zero lattice spacing (when the Dirac operator is antihermitian) correspond 
to exactly real eigenvalues at nonzero lattice spacing. Additional pairs of real modes appear 
for increasing lattice spacing. However, the number of spectral flow lines, \\{m), with an 
odd number of real zeros remains the same. We will identify the number of such flow lines as 
the index of the Dirac operator and study Dirac spectra for a fixed index. In the continuum 
limit this index, by the Atiyah-Singer index theorem becomes equal to the topological charge 
of the gauge field configurations. In Random Matrix Theory, this index is determined by 
the block structure of the Dirac matrix. All results in this paper are derived for fixed index 
of the Dirac operator. It is of course also possible to sum over all sectors with a given index. 

In this paper we elaborate on and provide detailed derivations of results announced in 



the Letter |20| and the proceedings [28|]. One of the simplifying features of that paper 
was that double-trace terms in the chiral Lagrangian were ignored (which can be justified 
based on large- Ai'c arguments 291]). Here we compute their contribution to the microscopic 
spectrum directly from the chiral Lagrangian. We also show how double-trace terms can be 
included in the Wilson chiral Random Matrix Theory. We study in detail the a — )■ and 
four-volume — ?■ oo limits of the analytical results. Furthermore, it is shown that the low- 
energy constants of Wilson chiral Perturbation Theory are constrained by a QCD inequality 
which follows from the Hermiticity properties of the Wilson Dirac operator. This constraint 
coincides precisely with the requirement of preservation of the Hermiticity properties of the 
Wilson Dirac operator. The constraints found are consistent with the existence of an Aoki 
phase. One important message from the calculations presented here is that the values of the 
low-energy constants of Wilson Chiral Perturbation Theory can be accurately determined 
by matching the predictions for the eigenvalue distributions to lattice data. Furthermore, 
our results describe analytically the eigenvalues that may cause numerical instabilities [sol 
when approaching the chiral limit at finite lattice spacing with Wilson fermions. 
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The analysis of the Letter 20] has been extended in various other directions. We have 
obtained the distribution of the chirahty over the real eigenvalues of the Wilson Dirac op- 
erator which is a lower bound for the distribution of the real eigenvalues. We also have 
obtained an upper bound which converges to the lower bound for small a. In the same limit 
the expressions for the microscopic spectral density of simplify and can be generalized to 
a nonzero number of flavors. We discuss the distribution of the tail states in the gap as well 
as a comparison with the scaling properties of such states as found in lattice simulations. 
We also perform a saddle point analysis of our analytical result and obtain a simple explicit 
expression for the gap of the Dirac spectrum. 

The paper is organized as follows. After a brief review of relevant properties of the 
Wilson Dirac operator given in section [TTl we define Wilson chiral Perturbation Theory at a 
fixed index of the Wilson Dirac Operator in section IIIII The derivation of the microscopic 
spectrum is given in section IIVI It is followed by a discussion of various limits in section |Vl 
A chiral Random Matrix Theory for Wilson fermions is shown to reproduce Wilson chiral 
Perturbation Theory in the microscopic limit in section IVll Finally, before concluding, we 
discuss the relation between the low energy constants of the Wilson chiral Perturbation 
Theory and the Hermiticity properties of the Wilson Dirac operator. Some technical details 
are referred to Appendix A, Appendix B and Appendix C 



II. EIGENVALUES OF THE WILSON DIRAC OPERATOR 

We start with a discussion of general properties of the Wilson Dirac operator. Some 



of these results have been discussed extensively in the literature (see, e.g., refs. 
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22 



24 



|26|, |27[), but are included here to make this paper self-contained. 

The Wilson-Dirac operator will be denoted hy D = Dw+m. The lowest order a correction 



was introduced by Wilson, 



Dw = ^7m(v,. + v;) - \av;v,. (1) 



It is written in terms of forward (V^) and backward (V*) covariant derivatives, and m is 
the quark mass. 

The Wilson Dirac operator, Dw, is not anti-Hermitian at non-zero lattice spacing. It 
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retains only 75-Hermiticity: 



(2) 



The eigenvalues, A^, of D]^ are then either real or they occur in complex conjugate pairs. 
The 75-Hermiticity of implies that the operator 



D5 = 1b{Dw + m) 



(3) 



is Hermitian. Since at non-zero lattice spacings a the axial symmetry is lost {{Dw-, 75} 7^ 0), 
the eigenvalues of -D5 do not occur in pairs of opposite sign. 

The eigenvalues, A^(m), of are nontrivial functions of m. If \\{m^ = 0, then 



D 



(4) 



so that the real eigenvalues of Dw can be obtained from the zeros of A|(m) with the corre- 
sponding eigenf unctions of and Dw beeing identical at m = rric- If we act with on 
the normalized eigenfunctions (pj = \j) of Dw we obtain 



(5) 



For m close to 



750i = {j\l5\j)(pj + S4>, 

with {6(f)\j) = and 50 ~ 0(m — rric). Therefore to 0{m — rric 



We thus find 



22| 



A5(m) = (m - mc)(j|75|j). 



(6) 



(7) 



(8) 



When evaluated at m = rUc, ie. at A^ = 0, the slope is thus exactly given by the chirality 

= 0175|J). (9) 



dm 



Furthermore, it can be shown 



22| that the chirality of the eigenfunctions 
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(10) 



vanishes for complex eigenvalues of Dw but is generally nonzero for real eigenvalues. 

To compute the spectral density ps of we introduce a source z that couples to ip'j^ip. 
The operator entering in the fermion determinant is then given by 



(11) 



and the resolvent of is defined by 



G{z) 



Tr 



+ z — ie 



(12) 



where \\ are the eigenvalues of D^. The density of eigenvalues of then follows from 

P5(A^ m; a) = (> ; 5{\l - \') ) = 3lm[G(-A^)],^o. (13) 



\ k I 



The resolvent ( fT2l) is a partially quenched chiral condensate. We will compute it by means 
of the graded technique, where the unphysical 'valence' determinant is canceled by an inverse 
determinant after differentiation with respect to the source 5^. In order to be able to write 
the inverse determinant as a bosonic integral, it is essential that an infinitesimal increment 
is added to a Hermitian operator 16 1. 

We also consider the chiral condensate corresponding to the regularization introduced in 
Eq. ([12]) 

1 



E(m) 



Tr- 



(14) 



Dw + m- ie75 

Care has to be taken in order to relate E(m) to the spectrum of Dy/. The discontinuity of 
S(m) across the real axis is given by 

1 1 



1 
1 

2Tii 



Tr 
Tr 



_ {Dw + m) - ie75 {Dw + m) + ^€75 

75 75 



(15) 



/y5{Dw + m) - ie 75(1)14/ + "i) + ze_ 

This expression can be rewritten in terms of eigenvalues and the normalized eigenfunctions 
I A;) of (recall that for A^, = the eigenfunctions of D^ merge with that belonging to a 



6 



real mode of Dw) 

e {k\^,\k) 



^xM = 1(11 

\ k 



5^5(A^(m)) (fc|75|fc) 

k I 



^ ^ \d\l/dm\ 

where are the eigenvalues of D^. Using Eq. this can be written as 

Px(A^) = ( E ^(^^ + ^'^) sig^((^l75|fc)) ) . (17) 



This shows that is the distribution of the chiralities over the real eigenvalues of Dyy. 
Integrating this expression over we obtain 

[ pAX^)dX^ = (j2 sign((A:|75|A:)) \ . (18) 

This is the average index of the Dirac operator. The index for a given gauge field configu- 
ration is defined by 

iy= J2 sign((A;|75|A;)). (19) 

Since the real modes of Dyy correspond to the vanishing eigenvalues of we also consider 
P5(A^ = 0,m). A similar calculation as above results in 

Because |(A;|75|A;)| < 1 we have the inequality (also valid for fixed index) 

Px(A'^) < Preai(A^) < p,{\' = 0,m = X"^ ; a) , (21) 

where 

f \W\ - 
PrcallA ) = 



E^(^r + A^)\. (22) 



As a special case the average number, (nrcai), of real eigenvalues of Dw for gauge field 
configurations with index v is bounded by 

y"rfA^p;;(A^) = z/ < Keai) < j dmpliX'^ = Q,m). (23) 
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III. WILSON CHIRAL PERTURBATION THEORY 



The operators contributing to Wilson chiral Perturbation Theory have been considered in 



a series of papers 
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4l3l]. It is constructed in terms of a double expansion in both the usual 



parameters of continuum chiral Perturbation Theory and the lattice spacing a. Depending 



on which counting scheme one uses, various terms contribute to a given order 3l|-l33|. We 
concentrate on the microscopic limit, the e- regime, in which the combinations 

mV, zV, aV (24) 

are kept fixed in the infinite- volume limit ^ — ?■ oo. To leading order in these quantities, 
the low energy partition function for lattice QCD with Wilson fermions then reduces to a 
unitary matrix integral. Up to the infinite- volume chiral condensate S and three low-energy 
constants determined by the discretization errors of Wilson lattice QCD, this integral is 
completely determined by symmetry arguments. 

Since we are computing quantities up to and including 0{a'^) effects, it is important that 
all effects to this order are contained in the on-shell effective Symanzik action including 0(a^) 



terms. This problem has been analyzed in detail by Sharpe for the p-counting-scheme 



19|. 



Here we reconsider the argument for the e-regime. One correction to relevant continuum 



operators corresponds to a'^ip{x)'j^'Vf^ip{x) [l9|. In the e-regime, such terms are suppressed 



by l/VV with respect to the leading terms and need to be considered only at higher orders 
in the expansion. Additional lattice artifacts of potential 0{a'^) in the form of contact terms 



have been analyzed .ef. M, and they are .nhlariy =upp.e.ed. 

As stressed by Leutwyler and Smilga [5|, the topological charge of gauge field configu- 
rations plays an important role in the microscopic spectrum of the continuum theory. The 
u zero modes of the Dirac operator distort the eigenvalue spectrum of the non-zero modes 
and lead to distinct predictions in sectors of fixed u. Likewise it is natural to study the 
microscopic spectral density of the Wilson Dirac operator for fixed index u. At the level of 
the chiral Lagrangian it is a priori far from obvious how to implement the projection onto 
sectors corresponding to a fixed index of Dw- Quite remarkably, such a projection can be 



achieved through a Fourier transform [20|: We decompose the partition function as 



Zn^ (m, z;a) =^ Z"^^. (m, z; a) (25) 
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with 



Z]; (m,2;a) = / dU defU e^^^^l (26) 

Ju(Nf) 



The Lagrangian C{U) is defined by 



C{U) = i(m + z)STr[/+ ^(m-z)ETrf/^ (27) 

Below we will demonstrate that in the microscopic domain this corresponds to an ensemble 
of gauge field configurations for which the index of D^/ as defined by Eq. (I19p is equal to v. 
The partitioning into sectors of fixed v corresponds at the level of the unitary group integral 
to writing an SU{Nf) integral as an infinite sum of U{Nf) integrals. 

In the Lagrangian (127|) . S is the usual infinite- volume chiral condensate while Wq, W-j 
and VFg are low-energy constants that quantitatively determine the discretization errors of 
Wilson fermions. The terms proportional to Wq and Wj are believed to be suppressed in 



the large limit 29| , and they were not included in ref . 20[ . Here we include the effect of 



these terms explicitly. To lighten the notation, we introduce the scaling variables 

m = mSV, z = zEV, a\ = a^WiV (28) 

with i = 6, 7, 8. 



IV. GRADED GENERATING FUNCTION 



To obtain the distribution, p^, of the chirality over the real eigenvalues of -Dvy and p\ at 
fixed index v we will use the graded method. The generating function in this case takes the 
form 

and we have two spectral resolvents 

Yy(vn\d) = - lim -^Z](r, , iii(m, m', 2; = 0, z' = 0; a), (30) 
m'^m am ' ' 

and 

G"'(z, m; a) = - lim , (m, m, z, z'; a). (31) 
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Note that the Nf physical quark flavors are not coupled to the source z. 

The presence of fermionic as well as bosonic quarks in this generating function gives rise 
to a graded structure ( "supersymmetry" ) of the corresponding chiral Lagrangian 

Z-^^^,^,{M,Z;a,) = JdUSdetiUy (32) 

^^i^STTiMlU-U-^])+i^STT{Z[U+U-'^])+dllST:t(U-U~'^)]^+^[STT(U+U 

where Ai = diag(m . . . m, m, rh') and Z = diag(0 . . .0, z, z'), with lm{z') < 0. The integra- 
tion is over the maximum Riemannian graded sub manifold of Gl{Nf + 1|1) 3 (see 0,13 
for notation and more on the graded method). For Nf = we will use the parametrization 

/ e*^ \ [ a\ 
U = { exp . (33) 

\^ y \f3 J 

If we could derive this partition function directly from the microscopic theory we would 
end up with convergent integrals. Relying only on symmetry arguments is not sufficient 
to obtain the partition function, and it can only be equivalent to the microscopic partition 
function if the integrations are convergent. For fermionic integrals convergence is not an 
issue, but bosonic integrals necessarily have to be convergent. 



In [s] the graded partition function was evaluated for the ant i- Hermit ian Dirac operator 
with positive mass and dj = and 2 = 0. The term relevant for the convergence of 
the integral was exp(— mcoshs) so that the generating function is defined for rh on the 
imaginary axis with a positive real infinitesimal increment, and can be used to calculate the 
spectral density localized on the imaginary axis. We could also have calculated the spectrum 
of by including the source term 275. This would have resulted in the bosonic integrand 
exp[— rncoshs — f sinhs]. The corresponding integral can be continued analytically to the 
entire real z-axis by shifting the s-integration by s — ?■ s ± 7rz/2 with the sign determined by 
the imaginary part of z. This transformation can be accomplished hj U iU for Im(i) < 
and by [/ — 7- —ill for lm{z) > 0. The expression fl32l) is thus valid for lm{z) < 0. 

Next we consider the partition function with df ^ 0. Let us first consider the case that 
= 0, and we consider the partition function as an analytical function of dg. The 
partition function ( l32l) is analytic in dg for Redg > 0, whereas the partition function of Ref. 
6|] can only be continued to Redg < 0. Depending on whether we consider the spectrum 
of or 75-Dvi/ the partition function can be analytically continued to a different part 
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of the complex ag-plane. Below we will argue that the physically interesting parameter 
domain is Re Sg > 0, which is accessible for the integration contours of Eq. fl32|) valid for 
the evaluation of the spectrum of -D5. Given this parametrization, the partition function is 
defined in only a subset of the complex and aj plane. For real constants the condition is 
that Sg — 67 — Sg > 0. Notice that in order to correspond to the QCD partition function, 
the limit — t- should be regular. 

We will now explicitly work out the generating function f p2|) in the quenched case where 
Nf = 0. In order to do so we will use the parametrization Eq. (!33|) which has a fiat measure 
{g]. In this representation the generating function with rh > and Im(i') < becomes 
(after a and /3 are integrated out) 

Zi|i(m,m i';ai) (34) 

ds — e'-*^"'*^'' exp[— m sin(6') — zm'sinh(s) + izcos{6) — zz'cosh(s) 
■00 J~n 27r 

+4a^(-z sin(^) + sinh(s))2 + 4a?(cos(^) - cosh(s))2 + 2a^(cos(2e) - cosh(2s))] 
X ( — — sin(6') + i— sinh(s) + cos{6) + cosh(s) 
-4{dl + d^j){sm\e) + smh\s)) + 2a^(cos(20) + cosh(2s) + e'^+' + e"*^"")). 

The s-integral is convergent for cig — fig — > 0. In this case it is easily checked numerically 
that Zi\i{m,m, z, z] d) = 1, as required by the definition of the generating function. We 
discuss the convergence criterion in detail in section IVIIBI below. 



A. The Microscopic distribution of the chirality over the real eigenvalues of D^K 

In this section we show that the parameter u in the chiral Lagrangian is the index of 
the Dirac operator Dw- Since we already have computed the generating function ( l32l) the 
condensate defined in Eq. (IT^ can be expressed as 

T.'^{m;di) = — lim -^ZJ^'h (m, m', z = 0, = — ze; Cj). (35) 

rh'-^m dm ' 

This corresponds to a resolvent regulated with ie'-f^ instead of ie, and as was discussed in 
section II, it has quite different properties with a discontinuity that does not give the spectral 
density but rather the distribution p'f of the chirality over the real eigenvalues of D\y given 
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in Eq. ([T7D. It follows from §^ that for lm{z' = z) < 

/CO PTT ^0 

ds — sin(6')e'-*^~*'''^ exp[— m sin(6') — zm sinh(s) — e cosh s 

+4a^(-i sin(^) + sinh(s))2 + 4a^j{cos{e) - cosh(s))2 + 2a^(cos(2^) - cosh(2s))] 

(771 777 
- — sm{e) + i— sinh(s) - 4(ag + d'^){sm'^{9) + sinh^(s)) 

+2dl (cos(2^) + cosh(2s) + e'^+' + e"*^"^) + -) . (36) 

The discontinuity of this resolvent is equal to microscopic limit of the distribution of the 
chirality over the real eigenvalues of defined in Eq. ([1] 



P^(A'^;a.) = (y1 5(Af + A^)sign[(A;|75|A;)]\ = -Im[S'^(m = A^; a,)]. (37) 



It can be verified numerically that with the resolvent (136|) 

p;^(A^;a,)dA^ = z/. (38) 



This demonstrates that the sectors introduced in Eq. (126|) correspond to a Dirac operator 
with index u, as defined in f|T9l) . 
In Fig. [1] we plot the distribution of the chiralities over the real eigenvalues of D^, 
p'^{X^ , di), for u = 1, 2, 3 for Sg = 0.2 and 66 = 07 = 0. For this value of ds additional pairs 
of real eigenvalues appear only rarely and p^(A^, di) is a good approximation to the density 
of the real modes of Dw- Indeed, the presence of u real eigenvalues is clearly visible. As we 
will show below, we can make an even more precise description of this in the limit of very 
small dg. In that limit we have u real modes which turn out to be described by the x 
Gaussian Unitary Ensemble of Random Matrix Theory. 



B. The Microscopic Spectrum of 

In the continuum (a = 0) the eigenvalues iXk of Dw and eigenvalues Xk of are related 
to each other through (A^)^ = (A|)^ — m^. In this case, the microscopic eigenvalue density 
of D,, 

pliX" > m,m;a, = 0) = J'l p'-(^(A5)2 _ ^2) (39) 
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FIG. 1: The quenched distribution of the chirahty over the real modes of Dw plotted for z/ = 1, 2 
and 3 with ag = 0.2. For the plot we have chosen ag = 07 = 0. The real modes repel each other, 
making the preferred locations of the eigenvalues clearly visible. 



follows from that of Dw- In the quenched case we have 



m 



P'W = ^ [JuW' - ^.+i(A)4_i(A)] . (40) 

To obtain the eigenvalue density, pg, of the Hermitian Wilson Dirac operator for 
non-zero values of a we first evaluate the resolvent 

G^{z,m;di) = — lini — Zj^M (m, m, 5' — ie; cij) (41) 

to find 

/oo PIT rn 

ds / — icos(^)e('^-^)" (42) 
■00 i-TT 2vr 

X exp[— m sin(6') — im sinh(s) + iz cos{6) — i{z — ie) cosh(s) 
+4a^(-z sm{e) + sinh(s))^ + 4a?(cos(^) - cosh(s))2 + 2a^(cos(2^) - cosh(2s))] 
X ( — — sin(^) + i— sinh(s) + cos{6) + cosh(s) 

-4{dl + d^^){sm\e) + smh\s)) + 2a2(cos(2e) + cosh(2s) + e*^+" + e"'^"^) + ^). 

The microscopic spectral density of -D5 = 75(1^14/ + m) in the quenched limit can then be 
expressed in terms of the imaginary part 

pl{X\m;di) = ^Im G''(-A^m;a^). (43) 
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FIG. 2: Left: The microscopic spectral density of for = 0, 1, 2 and 3 with, m = 3, ag = 0.2 
and ag = a? = 0. At nonzero value of the lattice spacing the zero modes spread out into a region 
around x = m. For negative values of v the spectral density is reflected at the origin. Right: 
After summation over z/ with a Gaussian weight and integrating from = up to = x. That 
is, the average number of eigenvalues below x. For the plot we have chosen (z/^) = 1. The values 



of m and are as in the left figure. Compare with Fig. 1 of 



In Fig. [2] this density is plotted for four values of u for fixed d and rh. The similarity of 
p^(A^) and p'^{X^) for A^ ~ A^ — m (see FigurelUfor a direct comparison) is not accidental: 
As we show in section |Vl for ag ^ 1 the z/ real modes of Dw are mapped directly to u modes 
in the vicinity of m, see Eq. flS^ . 

It is of course also possible to sum over all sectors with a given index, for an example see 
the right hand panel of [2J 



In Section II we discussed the mass dependence of p'^{X^ = 0,m; Sj). The density of the 
real eigenvalues of Dw satisfies the inequality fl2T]) . In Fig. [3] we plot both sides of this 
inequality function of A . 



C. The effect of Wq and Wj 



In the derivation above we have explicitly included the effect of Wq and Wj in the La- 
grangian and then performed the relevant super-traces and fermionic integrations. Here we 
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p^''(A^;a) (solid curve) 
P5^(X^=0,m = X^;a) (dashed curve) 
a = O.Z (black) 
a = 0.4 (red) 




FIG. 3: The distribution of the chirahty over the real eigenvalues of Dw p'^^{X^; en) (solid curves) 
and P5^'^(A^ = 0, m = A'^; en) as a function of A^. The values of ag (here denoted by a) are given 
in the legend of the Figure (ag = 07 = here). The density of the real eigenvalues of D\y is 
bounded between the solid and the dashed curve. 

point out an alternative way to compute the effect of and W-j. Not only will this give us 
a direct way to quantify the effect of and W-j on the spectrum of the Dirac operator, it 
also provides a simple way to include the effects of these terms in the chiral Random Matrix 
Theory discussed below. 

Since Wq and Wj are coupling constants of double-trace terms, we can re-express the 
microscopic partition function (l26l) as 

„ 1 



2 2 
«6 y? 



dyedyr e ^^i^ei ^^i^^i Z]^^,(m - ye, z - yr, 0, 0, ag)- (44) 



Here we have written the expression valid for negative values Wq and Wt. If Wq > the 
shift of rh is instead along the imaginary axis, rh — iy^. For Wj > we analogously shift z 
along the imaginary axis, z — iyj. 

Since exactly the same rewriting is valid for the generating function fl32|) also the quenched 

spectral density at non-zero values of Wq and Wj follows from that with Wg alone 

1 /"oo _ fi _ y? 

dyedyje ^ ^ p^(A5 - ^7, m - yg; 0, 0, ag). (45) 



P5(A ,m;a6,a7,a8) 



IGiraQdj 

Interestingly, for Wq < the quenched spectral density of D^j is a Gaussian average of 
spectra with a smooth distribution of quark masses. The integrals over y^ and yj in 
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and (H5|) can only be interchanged with the noncompact integral in the partition function if 
the inequality cig — Sg — > is satisfied. Furthermore, the shift must be along the real axis 
in order that the discontinuity across the real axis remains linked to the eigenvalue density. 

V. LIMITING CASES 

In this section we discuss various limiting cases of the spectral density of and Dw- 
Since the Dirac spectrum at W^j 7^ is given by a Gaussian integral over the Dirac spectrum 
for Wq = Wj = 0, only the dependence on Wg will be analyzed in this section. 

A. Small dg-limit for z — rh fixed 

We first show explicitly that in the limit Sg — the microscopic resolvent ( H2l) reduces 
to the known analytical result of the continuum limit. To derive the cig — t- limit of the 
microscopic spectral density, we use the integrals 




rh sin 8+iz cos 6 



( 



z + m 



z — m 




sign(z)A/^^ — Tn- 



(46) 




(47) 



valid at fixed \m — z 



7^ 0, and obtain 



lim G"^(i, m; as) 



-\K,{vm + m)(/32/,+2(r/) + + - m)(/3"'/.-2(^) + IM)] 

~[l3h+M + l^~^h^MW + rri)(3K,.i{T^) + (£ - m)/3-^JC+i(r/)] 

+ '-KMW.+M + p-'h-Ml (48) 



with 



7] 





(49) 
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We now use the recursion relations 

V 

K,+i{v) = K,_,{r]) + ^Kiv), (50) 
and the Wronskian identity to obtain the resolvent of in the limit ag — ?■ for z — rh fixed 

G^{z,m) = -z{h{r])KM + h+^K.^M) + ^r^- (51) 

z — m 

As we will now demonstrate, this is the resolvent of for a = 0. The resolvent of Dw at 
a = can be expressed as 

The eigenvalues of are X\ = ±\/ (A^)^ + for ag = resulting in the resolvent 



D5 + z / \ ^ — (AF)^ — / z — m 
(Ar)2 + m2 -^2 / + 



-2 



z — m 

x=i\' z'^—m? 



(53) 



Inserting the Sg = result for S'^, see eg. Eq. (45) of |6[, we consistently find ( ISTl) above. 
We conclude that in the limit ag — )■ the eigenvalue density p5(A^,m;ag) is given by the 
standard Sg = result f l39|) as long as |(A^ — m)|/ag ^ 1. 

B. Small ag-limit for {z — rh)/as fixed 

As we have just seen, for a = the contribution of the zero modes of Dw to the resolvent 
of Z^5 is given by 

(54) 

z — m 
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For a 7^ the zero modes are smeared out so that the singularity in the resolvent is 
smoothened, and we expect the imaginary part of the resolvent to behave as 

ImlCiz, m)] r^l F{z-m), (55) 
as 



with F a peaked function. In Appendix A we extract the analytical expression for F from 
the small as limit of the microscopic result for the resolvent given in Eq. ( l42l) . Amazingly, 
the result is given by the eigenvalue density of a |z/| x |z/| Hermitian Random Matrix Theory 

„/,c; , 1 , sr^ Hi.(u)Hk(u) , . 

p^(A^m;as«l) = — exp(-n^)^ ^^^^^ ^ (56) 



1 e-"' 



4a8 1-{y-\)\^ 

with 



[i/,'(n)-if,+i(u)if,_i(n)], 



u = (57) 



This is the familiar spectral density of the u x u Gaussian Unitary Ensemble shifted by m 
and rescaled by l/4as. For u = 1 the result is a simple Gaussian. 

For the normalization of the Hermite polynomials we used the convention 

/oo 
dx Hj{x)Hk{x)e-''' = Sjk'i^kly/^. (58) 
-oo 

The l/ag contribution to the spectral density is therefore normalized to v. For z/ = the 



leading small-ag corrections are 0(log(a8)), see [Appendix A| for details, and will not be 
considered to the order we are working. 



C. For small ag the distribution is given by 

The distribution of the chirality over the real eigenvalues of Dw, p^(A^; Og), is obtained 
from the resolvent of -D5, Eq. (l42l) . by replacing the pre-exponential factor icosO — — sin^ 
and putting z ^ 0. As discussed in section IIV Al to leading order in dg, this is equal the 
spectral density of the real eigenvalues of Dw- For small dg the width of the spectrum ~ dg 
and we can thus consider the limit of small dg with m/dg fixed. This is exactly the limit 
that was considered in previous section. 
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To leading order, only the negative exponent of — sin 9 has to be taken into account 
(which, up to a minus sign, is the same as the negative exponent of icosO). We thus find 



Preai('^ ~ ''T^'-i (^s) = P^i^ ,'fh]a^) for ag — 7- and (m — A )/as fixed. (59) 

This fact is also demonstrated graphically in Figure |H The two distributions merge because 
the real modes are almost chiral. To see this we start with a real mode of Dy/ 22 1 

Dw4>j = Af (60) 

It follows that 

^50i = 75(^14^ + m)(l)j = {Xf + m)'y5(j)j. (61) 

Now if the real modes of Dyy are chiral 750^ = ±0j then the u real eigenvalues of Dyy are 
mapped onto u eigenvalues of with a trivial shift by m. More precisely, 

750i = (j|75|j)0i + with = and 50 ~ 0(a), (62) 

so that 

^50, = ira + Af )(j|75|j)0, + 0(a). (63) 

Since the distributions of the two merge in the small ag limit, (see Eq. ( l59i) ). this explicitly 
confirms that the chirality of the real modes is unity to leading order in dg. 



D. Scaling of Smallest Eigenvalue 

For 1^ = 1 and ag ^ 1 the distribution of the single real eigenvalue of Dw takes the 
Gaussian form 



Pr.:;(A»';a,«i) = j-^exp^^-^). (64) 



In physical units the width of the distribution is 



a = (65) 
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- p ^(A ,m;a) 
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FIG. 4: The microscopic spectral density of D5 for u = 3 with m = 3 and as = 0.1 and oe = a? = 0. 
Also plotted is the distribution of the chirality over the real modes of D\y shifted by the mass, that 
is /0^(A^ — m; as). Since the real modes are almost chiral for small as the two distributions merge. 



This is also the width parameter of the Gaussian tail of the spectral density of inside 
the gap. In lattice simulations of the Wilson Dirac operator, a scaling of the width of 
the distributions of the smallest eigenvalue ~ a/y/V has been reported 30| for Nf = 2 
simulations. 

It is also instructive to consider the ratio of the width parameter and the average spacing 
AA = 71 /T,V of the Dirac eigenvalues at a = 0, 



— = (8a^iys^)--- = as . 

AA TT TT 

This gives an intuitive interpretation of the dimensionless low-energy constant 



(66) 



E. Mean Field Limit 

For large m, z, and the graded generating function can be evaluated by a saddle 
point approximation. This limit corresponds to the lowest non-trivial order in the usual 
perturbative expansion (j)- regime) as considered in [l^ . We will focus, in particular, on the 
behavior of the spectral density ps near the edge of the spectrum. At mean field level the 
dependence on the index of the Dirac operator is suppressed and we will therefore start from 
the u = expression and drop the index u below. Since the Dirac spectrum at Wqj 7^ is 
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given by a Gaussian integral over the Dirac spectrum for Wq = Wj = 0, only the dependence 
on Ws will be taken into account in this section. 

The expression for the spectral density can be written as 

-1 poo PTT jn J 

pr'{m,z;ds) = -Im / ds / e^/(™) e^^(-.^'^^B)p(m, f , ag) |,,_,, (67) 
vr i-TT 2vr dz' '^-^ 

where Sf, Sb and P can be read off from Eqs. f l4T|) . f H2|) . By shifting integration contours 
according to 

Q = ^(_r + -), 

TTZ 

s = t--, (68) 
the fermionic and bosonic exponents become real and equal up to a sign 

Sf{r) = mcoshr + zsinhr — 2ag cosh2r, 

Sb{t) = — mcosht — zsinht + 2ag cosh2t. (69) 

Also the prefactor in terms of these variables 

P{s,r) = — [mcoshr + zsinhr + rncosht + zsinht — 2ag((e^ + e*)^ + (e~^ + e~*)^)] (70) 

becomes real (an overall factor of i is included in the fermionic integration). 

For large m and z the integrals can be evaluated by a saddle point approximation. It is 
convenient to introduce u = sinh s as new variable so that the potential (!69|) is given by 



Sb{u) = -mvTT^ - zu + 2dl{2u^ + 1). (71) 

The edge of the spectrum is the point where two real saddle points coalesce and move into 
the complex plane. At this point 

Si{ug)=0, and 5^'(%) = 0, (72) 

where Sb{u) ia given in Eq. fITT]) . The second equation results into 

2/3 



m 



which combined with the vanishing first derivative leads to the mean field result for the gap 
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The mean field spectrum is symmetric around the origin and above we displayed the positive 
solution for the gap. 

Depending on the position of with respect to the position of the gap Zg we can distin- 
guish three parameter domains: 

i) |A^| < Zg, then all saddle points in terms of the r and t variables are real. The leading 
saddle point determines the fermionic integral but does not contribute to the imaginary 
part of the bosonic integral (which gives the imaginary part of the resolvent). The 
imaginary part is given by a subleading saddle point when combined with the fermionic 
contribution gives an exponentially suppressed tail. Combining this with the fermionic 
integral gives the spectral density. In the small ag limit, when coshs ~ sinhs at the 
saddle point, the bosonic integral becomes Gaussian resulting in the Gaussian tail 

^-2(m-8al)-^6a^J- ^^^^ 

ii) |A^| > Zg, then a pair of real saddle points has turned into a pair of complex conjugate 
saddle points. Only one of the saddle points can be accessed by the bosonic integration 
contour, but they both contribute to the fermionic integral. When the saddle points 
of the bosonic and fermionic integrals are the same, the exponents cancel resulting in 
a smooth contribution to the spectral density. When they are different, they result in 
an oscillatory exponent which is subleading because of the prefactor. It is instructive 



to work out the case ag = which is done in Appendix B Also the non-zero ag case 
is discussed in this appendix. 

iii) A^ ~ Zg, then two saddle points are close and the exponents can be approximated by 
a cubic polynomial. This is the scaling domain where V{z — ZgY^"^ is kept fixed in the 
thermodynamic limit for fixed m and a. In principle, the exact generating function 
can be evaluated in this limit, and according to universality arguments it should give 
a spectral density that can be expressed in terms of Airy functions as (see Eq. f ll45p 
for the definition of A) 

P5(x) = ^(Ai'(x)2 - xA:i^{x)), with x = {zg - z)/A. (76) 

The leading order asymptotic behavior of the Airy function on both sides of Zg follows 
from a saddle point approximation of the generating function in this domain. This is 



shown in [Appendix B 
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■ Airy soft edge 

Exact microscopic 




FIG. 5: The Airy behavior of (A^, z) dots for m = 30, = 1 and close to Zg vs. Eq. (1 76 p . 



F. Edge Scaling 

As discussed in previous section, the supersymmetric generating function can be evaluated 
in a counting scheme where m ^ 1, Sg ^ 1 and {z—ZgY^'^ ^ 1 are all of the same order. This 
is a universal scaling domain where the average spectral density is given by the expression 
in terms of Airy functions (see Eq. f l76|) ). We illustrate this in Fig. |5] for m = 30 and Sg = 1. 

The distribution, Pmm, of the smallest positive eigenvalue for the corresponding density 



Eq. ( 176|) is also known. It is given by the Tracy- Widom distribution 36 



31 for /3 = 2, 



Pn.in{z) = -^F2{{Z-Zg)/A), (77) 

F2(x) = e-^--^"-^)^'^^)'^^ (78) 

where q{x) is the solution of 

q"(x) = -xq{x) + 2q^{x) (79) 
with boundary condition that q{x) — > Ai(— x) for x — )■ —oo. 



VI. RANDOM MATRIX THEORY FOR THE WILSON DIRAC OPERATOR 



A chiral Random Matrix Theory for lattice QCD with Wilson fermions is constructed 
from the most general 75-Hermitian matrix. This random matrix Wilson Dirac operator has 
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the block structure 

( aA iW ^ 

Dw= \ ^ , (80) 

yiW^ aB 

where 

A = A^ and = B (81) 

are (n + z/) x (n + z/) and nxn complex matrices, respectively, and W is an arbitrary complex 
{n + v)xn matrix. In Eq. (|80|) and below we use tildes to indicate quantities in the Random 
Matrix Theory which are analogues of those in the quantum field theory. We relate the two 
sets in Eq. (1551). 

We take the matrix elements to be distributed with Gaussian weight 

V{A, B, W) ^ e-f Tr[A2+B^]-f Tr[14^iyt]^ (82) 

where N = 2n + u. Because of universality, results in the microscopic domain should not 
depend on the details of the probability distribution j4|, but for simplicity we take a Gaussian 
distribution. The partition function of the Wilson chiral Random Matrix Theory is then 
defined as 

Nf 

Z^^^ = I dAdBdW JJ dei{iDw + + 5^75) V{A, B, W). (83) 
^ /=i 

The matrix integrals are over the complex Haar measure and 75 = diag(l, ••• ,1,— I,-- - ,— 1) 
is a diagonal matrix with n + z/ diagonal matrix elements equal to 1 and n diagonal matrix 
elements equal to —1. 

For large quark mass the matrix = 75 (D^ + m) has n + u positive eigenvalues and 
n negative eigenvalues, whereas for large negative mass -D5 has n + u negative eigenvalues 
and n positive eigenvalues. Therefore at least |z/| spectral fiow curves of -D5 have to cross 
zero at least once. Since each crossing corresponds to a real eigenvalue of Dw we conclude 
that has at least |z/| real eigenvalues. The block structure of the matrix (180|) guarantees 
that the Random Matrix Dirac operator has index u. 

A chiral Random Matrix Theory for staggered fermions at finite lattice spacings was 
introduced in jsg]. 
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FIG. 6: Scatter plot of the eigenvalues of Dyy with i/ = 5 as obtained from Wilson chiral Random 
Matrix Theory. The real eigenmodes (black dots) are clearly visible. The width of the strip of 
complex eigenvalues in the thermodynamic limit equals |Re(z)| < Sa^Wg/Ti, as also follows from 
the chiral Lagrangian. See Eq. (j88p for the match between the parametrs of Wilson chiral Random 
Matrix Theory and those of the chiral Lagrangian. 

A. Prom Wilson chiral Random Matrix Theory to Wilson chiral Perturbation 
Theory 

We now consider the microscopic limit of the chiral Random Matrix Theory for Wilson 
fermions in which 

m~Ar-\ 5~A^-\ ar^N-^^^. (84) 

We have studied numerically the quenched eigenvalue spectrum of the random matrix Wilson 
Dirac operator, D^r, in this scaling regime. We find that the number of real eigenvalues of 
Dw is equal to i/ + 0(a^) for u ^ and z/ + 0(a^) for z/ = 0, while the remaining eigenvalues 
come in complex conjugate pairs, see Fig. |6l Moreover, we have found perfect agreement 
with the analytical predictions for in Eq. ( P7|) and in Eq. this is illustrated by 
Fig. [71 This agreement strongly suggest that the chiral Random Matrix Theory for Wilson 
fermions introduced above is in the same spectral universality class as the leading term of 
Wilson chiral Perturbation Theory in the e-counting scheme. 

To demonstrate this universality we will now show that the Wilson chiral Random Matrix 
Theory partition function reduces to that of Wilson chiral Perturbation Theory to leading 
order in the e-counting scheme. The explicit calculation will allow us to identify the low 
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FIG. 7: The average spectral density of the hermitian Random Matrix Dirac operator with 
n = 50 and 1^ = 1 over an ensemble of 100,000 matrices (marked by crosses) and the prediction 
from Wilson Chiral Perturbation Theory. 

energy constant Wg in terms of the parameters of the Random Matrix Theory. Inclusion in 
the Random Matrix Theory of the two other terms in the chiral Lagrangian to this order that 
couple with Wq and Wi will be discussed in the next subsection. For notational simplicity 
we keep the quark flavors degenerate. 

The first step is to express the determinants in Eq. fl83l) as Grassmann integrals. Then the 

nn 

average over the matrix elements of W can be performed by completing squares, see |3| 
for details. The four-fermion terms can be decoupled by means of a Hubbard-Stratonovitch 
transformation, leading to the partition function 



Z"^ (m,z;a) 



dQidQ2dTdT^ 



(85) 



X 



det^/2+^(mgi + T + m + z)det''/\t~aQ2 + T^ + m- 5)e-f ^^^^^-f Tr[Q?+Qi]. 



Here, Qi and Q2 are Hermitian Nf x Nf matrices and T is an arbitrary complex Nf x A^^^ 
matrix. 

Up to now this is an exact reformulation of the partition function ( l83l) . In the limit of 
large matrices valid for the scaling f l84p the integrals can be evaluated by a saddle point 
approximation. The saddle point manifold is given by 



T = UV~\ 



(86) 
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After absorbing V ^ in U and expanding the exponent to order m, z and we obtain 

Z"^ (m,5;a) = / dU [ dQidQ^ det"(f/) gf (rn+5)TrC/+f {m-i)wt 

mTr(Qi(7+Q2(7-i)+f a2Tr[QiC/-iQiC/-i+Q2[/Q2C/]-f Tr(Q2+Q2) 

The integrals over Qi and Q2 are Gaussian and are locahzed on the saddle point. The term 
of 0{a?) in the exponent results in higher order contributions and they can be omitted here. 
We thus obtain the partition function 

Z]^ fm, 5; a) = / dU det\U) (-+i)w+f (,n-i)wt-iv^TV(f/^+c/t^). ^g^^ 

JU{Nf) 

This shows that the Wilson chiral Random Matrix Theory partition function (!85|) reproduces 
the three leading terms in the e-expansion. The corresponding low energy constants are 
identified as 

Nm = mEV, Nz = zEV, ^ = a^WsV. (88) 
As discussed in section IVII CI this identification requires that Ws > 0. 



B. Double- Trace terms 



The double-trace terms of Wilson chiral Perturbation Theory are not generated by the 
Wilson Random Matrix Theory defined in Eq. (!83|) . Although such terms have been argued 
to be suppressed in large- A^^^ counting 29|, we would nevertheless like to be able to include 
them. As we will now show, this can easily be done. 

Let the Random Matrix Theory partition function be extended to 

/ , , "TTTTTT Tf5|-2 



ZUm, z; cie, dj, ag) = / dyQdy^ e ^^'^si Z'^. im -ye,z- yr, cig), (89) 

where the partition function inside the integral is given in (183|) . It follows from the discussion 
of subsection IIV CI that then also the trace-squared terms of Wilson chiral Perturbation 
Theory are reproduced by the Random Matrix Theory in the scaling limit (1841) . The partition 
functions flS^ lead to negative values of and Wj. Positive values of these constants can 
be obtained from fluctuations of the mass in the imaginary direction. 
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VII. HERMITICITY PROPERTIES AND THE SIGN OF Ws 



The sign and magnitude of the low energy constants We, Wj and Wg are essential for 
numerical simulations of lattice QCD with Wilson fermions. In particular the sign of Ws has 



been debated in the literature 



19 



31 



40|. The sign is important for understanding whether 



or not 
parity 



attice QCD with Wilson fermions enters an Aoki phase with spontaneously broken 
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42]. In this section we provide several arguments for why Hermiticity properties 



put constraints on the these low-energy parameters. To simplify the discussion, we restrict 
ourselves to the case where both Wq and vanish. 



A. A Wilson lattice QCD inequality and the sign of Wg 

With two degenerate flavors the determinant of the Wilson Dirac operator is positive def- 
inite and it is possible to formulate rigorous QCD inequalities. Based on such an inequality 
we argue here that for Wq = W-j = one must have Wg > in the Nf = 2 theory. Since the 
Wilson fermion determinant in lattice QCD is real, the measure with two degenerate flavors 
is real and positive 

det"^ {Dw + m) > 0. (90) 

As this is true for any gauge fleld conflguration, independent of the number of real eigenvalues 
of Dw, it follows that the two-flavor partition function, of lattice QCD with Wilson 

fermions in a sector with flxed number of real modes of is real and positive for all real 
values of m 

^]^;=2K^ = 0;a) > 0. (91) 

The overall sign of the partition function can of course be changed by introducing a multi- 
plicative constant. What the inequality states is that the partition function cannot change 
sign as a function of a real valued quark mass. 

The partition function in the e-regime must necessarily satisfy the same positivity bound. 
We note, however, that its sign depends on the index u of the chiral Lagrangian: 

^];,=2(^, ^; -al) = {-lyZ^^^^S^, im; al). (92) 
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FIG. 8: The two flavor microscopic partition function for a| = 1 (Wg > 0) and a| = — 1 (Wg < 0) 
plotted as a function of the mass for z = 0. Only the partition function with VFg > satisfies the 
QCD inequality. 



For odd u the partition functions with positive and negative have opposite signs at 
m = z = and, since for sufficiently large rh both partition functions have the same sign, it 
is not possible that they both satisfy the inequahty for real rh and z. As the plot in Fig. [8] 
demonstrates, the sign of Ws dictated by the inequality is positive. 

This can also can be obtained analytically from the chiral Lagrangian. The partition 
function with degenerate masses fl27|) can be expressed in terms of the eigenvalues of U. 
This allows us to rewrite the mass-degenerate partition function in the simplified form 



Z^^{m,z]d8) = det[Z''/^l^^{m,z;ds)]ij=i...Nf 



(93) 



where the one fiavor partition function is 

de 



27r 



^Qv ^pi cos((9)+jz sin(S)— 2a| cos(26) 



(94) 



In Appendix C| we show that for odd values of v the two-fiavor partition function with 
PVg > (lig < 0) is positive (negative) for m = and that for both signs of the 
curvature at zero mass is positive. 

Below we discuss the sign of W% from the perspective of the graded generating func- 
tion and from Wilson chiral Random Matrix Theory. In both instances we find > 0, 
independent of the number of fiavors. 
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B. Convergence and 75-Hermiticity 



The issue of the sign of the constants of Wilson chiral Perturbation Theory manifests itself 
also in the convergence properties of the graded generating functional. For Wq = Wj = 
the integrals in (I32p are divergent if Wg < 0. However, the alternative graded generating 
function 

Zr|i(A^, Z; al<0)= jdU Sdet([/)^ ^iSTr(A^[l/H-f/-]) + lSTr(2[C/-C/-])-aiSTr(C/^+C/-) ^95) 

is now convergent. This integral seems to provide a bona fide generating function for Wilson 
fermions with Ws < 0. In factjjt agrees precisely with the generating function suggested for 
regime calculations in ref. 19|. Nevertheless, the sign-issue has not disappeared, it has 
only resurfaced in disguise. To see this, we can compute the resolvent for the real eigenvalues 
of Dw with the alternative graded generating function 



— Z 

dm 



Z'(JM,Z = 0;dl<0) 



(96) 



By comparison with Eqs. (!32|) and (!95|) . it follows that 



d 



Z'(,,iM,Z = 0;dl<0) 



-I = 0,Z;a^>0) 



(97) 



z'=z=im 



dm 

Since the rhs is continuous when m crosses the real axis we conclude that for < the 
density of real eigenvalues vanishes. This is because the graded partition function, (^5]) . 
valid for < 0, is the generating functional of an ant i- Hermit ian operator with a spectral 
density on the imaginary axis. This would be an acceptable eigenvalue density if Dyy was 
anti-Hermitian, but this is only the case for a = 0. Away from a = the operator is only 75- 
Hermitian. So from the convergence of the graded generating function and the requirement 
of 75-Hermiticity of D]y it follows that Wg > 0. 

The possibility of a link between the integration domain of the graded generating func- 
tional in chiral Perturbation Theory and the Hermiticity properties of the Dirac operator 
is not new. In continuum QCD at non-zero quark chemical potential the Dirac operator is 
non-Hermitian for real values of the chemical potential and anti-Hermitian for imaginary 
values. Correspondingly, the graded generatin g fu nctional has two branches which differ by 



a transformation of the super Goldstone field 



44| . One is convergent for real values of 
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the chemical pot ential while the other is convergent for imaginary values of the chemical 
potential Q, 3- 

C. The sign of Ws and Wilson chiral Random Matrix Theory 

In this section we discuss the constraints on the constants in Wilson chiral Perturbation 
Theory due to 75-Hermiticity from the point of view of chiral Random Matrix Theory. 

In Eq. f l88p the coefficients of Wilson chiral Perturbation Theory where expressed in terms 
of the parameters of the chiral Random Matrix Theory, and it was found that Ws > 0. A 
chiral Random Matrix Theory where the sign of Wg is negative can be obtained by including 
an additional imaginary unit i in front of the chiral symmetry violating pieces of Dw in (|80|) , 
that is 

, iaA iW , , , 

Dws<o = \ , • (98) 



2 



iaB 



With this change, the Random Matrix Wilson Dirac operator, however, becomes anti- 
Hermitian rather than 75-Hermitian. We stress that this observation is in perfect agreement 
with the conclusion obtained from the convergence requirement of the graded generating 
function above: the spectral density of the anti-Hermitian Random Matrix ( I98p matches 
that generated by the graded partition function (!95|) which is convergent for Wg < 0. This 
has been verified numerically to high accuracy. Because Dws<o is anti-Hermitian, the deter- 
minant det{Dwa<o + has a definite sign for imaginary m in agreement with the chiral 
Lagrangian for Wg < and imaginary m. 

From the point of view of Random Matrix Theory, the different universality classes cor- 



47|. In the 



responding to whether one imposes 75-Hermiticity or not are well understood 
Random Matrix Theory literature this is referred to as Q-symmetry. 

Clearly the conclusions of this section have the same origin, namely that the case Ws < 
(and Wq = W-j = 0) is in conflict with the 75-Hermiticity of D]^ and the Hermiticity of D^. 
A positive sign of Ws is consistent with the existence of an Aoki phase. 
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VIII. CONCLUSIONS 



In this paper we analyzed the spectral properties of the Wilson Dirac operator in the 
microscopic scaling regime. Eigenvalues of the Wilson Dirac operator in this regime are 
responsible for spontaneous chiral symmetry breaking. Another feature of Wilson fermions 
is the existence of the Aoki phase with spontaneous breaking of parity. In the microscopic 
scaling regime one has analytical means for studying both. We have shown that the chiral 
Lagrangian for Wilson Chiral Perturbation Theory can be used to compute spectral prop- 
erties of the Wilson Dirac operator Dy/ and its Hermitian counterpart -D5. Using a graded 
extension of the chiral Lagrangian up to and including O{o?) effects, we have computed 
the microscopic spectrum of the Hermitian Wilson Dirac operator in the quenched the- 
ory. We have shown how to incorporate all possible terms of the Lagrangian to this order, 
including double trace terms. 

An alternative path to the microscopic spectrum goes through a chiral Random Matrix 
Theory based on symmetries of the Wilson Dirac operator. Here we have formulated such a 
theory and proved that the partition function coincides with that of Wilson chiral Perturba- 
tion Theory to leading order in the e counting scheme. Numerically, we have demonstrated 
that also the spectral correlation functions of the chiral Random Matrix Theory agree beau- 
tifully with the analytical predictions of Wilson chiral Perturbation Theory. 

Interestingly, we find restrictions on the possible values of the low-energy constants of 
Wilson Chiral Perturbation Theory arising from the imposition of 75 Hermiticity. It would 
be most interesting to have a more detailed understanding of this phenomenon, which seems 
to run deep. It is noteworthy that the bound we get from the chiral Lagrangian matches 
exactly with the bound we get, by an entirely different route, from our Wilson chiral Random 
Matrix Theory. In its simplest form, where we ignore double trace terms, it was found that 
this bound is consistent with the existence of an Aoki phase. 

The results presented here are are for the microscopic limit in sectors with a fixed index of 
the Wilson Dirac operator. The real modes of the Wilson Dirac operator have an interesting 
dynamics on their own. Remarkably, we find that very close to the continuum, the \v\ real 
modes behave like the eigenvalues of a \ v\ x \ v\ Hermitian Random Matrix Theory of exactly 
Gaussian weight. We emphasize that this is a universal result that has been derived from a 
chiral Lagrangian. In lattice gauge theory simulations it would be most interesting to study 
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just this subset of real modes of the Wilson Dirac operator. In particular, as we have shown, 
their dynamics carries detailed information about the low-energy constants of Wilson Chiral 
Perturbation Theory. 

In general, we suggest to use spectral properties of the Wilson Dirac operator to determine 
the new parameters of Wilson Chiral Perturbation Theory. While these parameters are 
unphysical, it is nevertheless crucial to establish their values if one wishes to extract physical 
observables from Wilson fermion simulations at finite lattice spacings. It will thus be of 
obvious interest to extend the first results presented here for the quenched theory to the 
corresponding theories with Nf light fermions. Work in this direction is presently underway 
51|. 

Acknowledgments: We would like to thank many participants of the Lattice 2010 Sympo- 
sium and the TH-Institute "Future directions in lattice gauge theory - LGTIO" at CERN for 
discussions. We thank C. Gattringer for clarifying correspondence. Two of us (GA and JV) 
would like to thank the Niels Bohr Foundation for support and the Niels Bohr Institute and 
the Niels Bohr International Academy for its warm hospitality. This work was supported 
by U.S. DOE Grant No. DE-FG-88ER40388 (JV) and the Danish Natural Science Research 
Council (KS). 



33 



Appendix A. SMALL ag-LIMIT OF RESOLVENT FOR {z - m)/as FIXED. 

In this appendix we consider the leading nontrivial expansion of for small ag at z/ > 0. 
We start from fH2]) with Wq = Wj = and then compute the small ag limit for [z — m)/a^ 
fixed. Our aim is to prove that in this limit the eigenvalue density is given by the familiar 
expression for the density of the Gaussian Unitary Ensemble with matrices of size u x h>. 

The fermionic ^-integrals in f l42p are all of the form 

= — [ rf^e^PV™'^^^+'^^°'^^+2aicos2e^ p = i/ - 3, . . . , i/ + 3 

27r 



— f (^^g*P^gv'm^-z^cos(6'+()!>)+2a| cos 29 
— r t/5)e»P(^-'/')e^"^^-^^™^(^)+2«8 COs2(e-<^)^ /gg\ 



with 



-i<f> 



1/2 

z — m \ iz m 

and cos cp = =, sm = :. (100) 



z + mj - ' yjm? - z 

To obtain the small ag and small m — z limit we expand the exponential functions in a Taylor 



series 



Zl = — r dee'p(^-*) y -}-(m^- z2)fc/2cos'=^(2a^cos2(e - </.))^ (101) 

The leading terms, which we will refer to as Z^, are obtained from exponents in the cosines 
with opposite sign to the sign of p 



2 _ ^2^|fc/2;^2^^2^^sign(p)0 

k+2l=\p\ 



y ^^(^2_^2^fc/2^2^ -^fcsign(p)0 



where we define 



k\l\1 

fc+2i=|p| 



i{z - sign(p)m) 

a = ^ . (103) 

2agy2 
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For p = we have Z^^q = 1. The sum above can be expressed in terms of Hermite 
polynomials, normalized according to (j58l) . for all p, 

H = (a8V^)H^^^i7H(a/(^v^)), (104) 

where the only dependence on the sign of p is in the argument a (noting that Hq = 1). 
We now turn to the bosonic s-integrals in ( 142|) . They are all of the form 

/oo 
ds e-P'e-'"'^^'^-''"'°''^^'^-'^'''^^ , for p = - 2, . . . , + 2 , 

■oo 

y-Pg-l (m+z)y- \ (z-m)/y-al {y^+y'^)^ (105) 



/o y 

after using the same definition as in fllOOp and changing variables y = e^. In the limit a — )■ 
with \rh — z\/ds fixed the leading contribution in an expansion in powers of dg is given as 
follows. 

For p > (p < 0) we change to rescaled new variables y = ds\/2/t {y = t/a^\/2) to 
obtain to leading order 

/"OO 

Zl = (agv^)"'*'' / dt tlPl-^exp[-at-tV2] 
Jo 



with Dp a parabolic cylinder function 48|, and a depending on sign(p) defined in f llOSp . 

For p = the behavior in cxg is different. Because of the exact result Z^^z = = m) = 
KQ{2a\) we expect a logarithmic singularity here. Splitting the integration in Eq. ( I105P into 
(0, 1) and (1, oo) and changing variables as for p > for the former integral, and as for p < 
for the latter, we arrive at 



^p=0 



dt 

;V2 i 



e 



ai(2aiA2+f2/(2ai)) 



(107) 



Because the saddle point is outside the integration contour the leading contribution comes 
from the lower endpoint of integration dg\/2, and we arrive at 

= - log(2a^) exp[-z(m + z)/2 - i{z - m)/2]. (108) 
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By combining the fermionic and bosonic integrals we obtain 

G'^(5,m,a8) = '-]^ztl^^im + z)[Zl^, + Zl]+'-iz-m)[Zl + ZU 

+al[Zl^s + zUi + zLi + Z^]} 
+{Zl^, + Zl_,) |^(m + z)Zt^, + '-{z - m)Zt_,, + al[ZU + ^^+2]} 
+2alZl_,[zU^ + Zl] + 2alZl^,[zU + Zl]] . (109) 

To obtain the leading small Og limit we simply replace the partition functions by those with 
a caron. 

For z/ > 3 all bosonic integrals have a positive subscript p > and no log-terms appear. 
For the fermionic integrals for > 3 it is immediately clear that in the limit ag — )■ only the 
fermionic partition functions with the lowest index have to be taken into account in each 
sum. The cases u = 0,1,2 have to be checked separately. Here special attention has to be 
paid to the fact that for sign(p) < we have a ~ 1/ag. 

On the fermionic side the term Zl_^^ + Zl_^i + Z^^^ + Zl_^ is sub leading both for z/ = 1 
and u = 2 whereas Z^ + Zl_2 is subleading for = 1. (This is also the case for z = 
relevant for the resolvent of the real eigenvalues of D^.) Below this is automatically taken 
care of by the factors {u — 1) and (z/ — 2). 

For u > the expression for the resolvent thus simplifies to leading order to 

G^>°(f,m;ag) = (110) 



a 



H^^i{-^)[D^^{a) + z/aL)_i_^(a) + u{u + l)D_^_2(a)] 



2 2(-i)/2ag2v^ L tV2 

_^ v^(z. -l) ^^ a ^ 2uD^^_,{a)] + ^^^^^4^^if._3(T^)/^_.(a) 

I zV2 i iV2 

For z^ = we obtain no order 1/ag terms, and the leading order is in this particular case 
given by 

i log(2a^)(z + m) exp[-z(m + z)/2]. (Ill) 
8 

Using the recursion relation 

Dp+i{z) - zDp{z) + pDp_i{z) = 0, (112) 
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one can simplify the first term in Eq. flllOp 

G''>%z,m;as) (113) 



zv2 



/—ex (X 
^V2 «v2 

For p a positive integer, the parabohc cyhnder functions can be written as 

D-M) = ^Jl ^J,y2^p[^y H,^ii^^ + P,.,{a)e-'/\ (114) 

where we define the following polynomials that have parity k 

lk/2] 
1=0 

with real coefficients ai. This relation follows from induction, using the recurrence relation 



f lll2p as well as the following 48 1 



D^{x) = 2-"/2e-^'/4if^^(x/v^) , for n = 0, 1, 2, . . . 
D^i{x) = ^e'''^eiic{x/V2) . (116) 

By inspection it is clear that the terms containing the polynomial Pk do not contribute to the 
imaginary part of the resolvent. Because of erfc(x) = 1— erf(x) and erf(a:) being odd, only 
the term proportional to unity contributes to the imaginary part. After using the recursion 
relation for the Hermite polynomials, 

Hn+i{x) = 2xHn{x) - 2nHn-i{x) , (117) 

we find for the imaginary part of the resolvent and hence the eigenvalue density of (with 
ag — and \z — m\/a^ fixed valid for > 0) 

pl>%z,m-a,) ^ -\m[G^>%z,rn;a,)] (118) 

TT 

-\Hl{a/{iV2)) - H,^M/{iV2))H,^r{a/{i^))]. 



y/^a^2^+^{v- 1)! 

This is the familiar density of the v^v Gaussian Unitary Ensemble shifted by rh and rescaled 
by l/4d8. For i/ = the density is simply zero at the same order. 
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Appendix B. MEAN FIELD LIMIT 



In this A ppen dix we give more details for the mean field results discussed in section IV El 



fsee also 



16 



49|). 



A. Mean Field Analysis of Graded Partition Function 

For large miV, zN and a^N , the integrals in the expression for the generating function 
of the resolvent can be evaluated by a saddle point approximation. Unless the Grassmann 
integrals vanish at the saddle point (we will see below that this indeed may happen), the 
leading order result can be obtained by putting the Grassmann variables equal to zero so 
that the integral factorizes into a compact and a noncompact integral which each can be 
approximated by a saddle point integral. 

A priori we can choose either the compact or the noncompact integral to derive the mean 
field result for the resolvent. This limit corresponds to the physical limit of taking the 
thermodynamic limit at fixed lattice spacing. 

The mean field approximation to the supersymmetric generating function is thus given 
by (provided that the Grassmann integrals do not vanish at the saddle point). 



Z(m, 5, z'] Sg) = Z^{m, z; a^)Z^{fh, z'; Sg) 



;ii9) 



resulting in the spectral density 

P5{z,m] as) = Im 



vr dz' 



Let us start with the fermionic integral given by 

de 



Z^ 



z, m: 



msin d-\-iz cos 6+2a^ cos(2S) 



TT 



(120) 



(121) 



For real z, the imaginary part of the integral vanishes so that the spectral density is given 
by 



P5{z,m; ds) = Z^{z,m;d8) --^lmZ^{z ,m; ds] 

71 az 



The fermionic integral can be rewritten as 



Zf{z,m; ds) = -Re f' ^^e-'^^+*^v^+2ai(i-2,2)_ 



(122) 



(123) 
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The extremum of the real part of the exponent is at ?/ = — m/SSg, so that for SSg < m, the 
integral can be approximated by expanding about y = —1. In the thermodynamic limit we 
arrive at 

Jo 

m-2al-P /2{r'h-So?) 

- — ■ (124) 

^Jm - Sag 

The ^-dependence in the exponent is subleading in the above expression, so that the 
fermionic part of the partition function does not contribute to the resolvent to leading 
order. For ag = the integral in Eq. ( I12ip is given by (see Eq. fj46l) ) 

y/ 777, ^ — 

Zf{z, m- ds = 0)= 2Io{V^f^^^) ~ ^ J"..^ (125) 
For z <^'m this expression is approximated by 

^rh—z^ /2rh 



in agreement with the asymptotic result (11241) for a = 0. 
Next we consider the bosonic integral given by 



(126) 



Z\z,m-a8)= ^5e-*msinh«-*i cosh .-2ai cosh 2s _ ^^3?) 

J —00 

The imaginary part of the partition function can be written as 

77J,' ci ) — — / l^girh sinh s+ii cosh s—2ag cosh 2s g—im sinh s—iz cosh s—2ag cosh 2sj (X28) 

2^ J —00 

It is convenient to shift the integration contour by —111/2 so that the exponent becomes real. 
The saddle point approximation to the bosonic integral is then given by 

^P(i, m; as) = (z^^) '"'^'"^^ (^29) 

where 

Ss{s) = —m cosh s — z sinh s + 2ag cosh 2s, (130) 

and the sum is over the saddle points. For z < Zg the saddle points are real and only 
the saddle points with S'^{sp) > contribute to Im(Z^). There is one saddle point with 
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S'J{sp) < which gives the real part of and is of the order of {Zf)~^ so that the graded 
partition function is normahzed correctly at the mean field level. We thus find that for 
z < Zg the real and imaginary parts of are determined by different saddle points. 

Ai z = Zg the real saddle point with negative curvature merges with the real saddle point 
with positive curvature (which determines the imaginary part of the partition function), and 
for z > Zg they turn into a pair of complex conjugate saddle points. However, only one of 
the two saddle points is accessible resulting in a bosonic partition function with a real and 
an imaginary part that are both determined by the same saddle point. 

Since the negative of the fermionic exponent is obtained from the bosonic exponent by 
replacing s — > iO, the saddle points for the bosonic and fermionic integral are the same but 
in the fermionic case both saddle points of the complex conjugate pair contribute to the 
partition function resulting in a real expression. 

Notice that replica symmetry or supersymmetry is broken for the contribution to the 
tail. The fermionic integral is always real and the imaginary part is due to the bosonic 
integral 50|]. The resolvent that can be derived form the fermionic partition function does 
not have an imaginary part for any number of flavors. The replica trick therefore fails for the 
fermionic partition function even at the mean field level. To get the correct result we have to 
select one of the two saddle points. This is the case for the bosonic partition function where 
only one of the two saddle points is accessible by deformation of the integration contour. 

A leading order saddle point approximation for the imaginary part of the bosonic partition 
function is accurate for a large parameter range. In particular, for large rh and z it covers 
both large and small Sg. There are two parameter domains where the derivation simplifies. 
First, for (jh — z)/a^ ^ 1, then coshs in Eq. f ll30p can be approximated by sinh(s). Second, 
for z close to the edge of the spectrum, where two real saddle points are close and the 
potential f ll30p can be approximated by a cubic potential. We first discuss the small ag 
limit. 



B. Tail of Eigenvalue Distribution for {rh — z)/a^ » 1. 

It is convenient to introduce u = sinh s as new variable so the potential (I130p is given by 
Eq. ( 17T|) . For {m — z)/as ^ 1, the leading saddle point 1, so that yl + v? ^ For 
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z > the leading saddle point is negative (and mutates mutandis for £ < 0), so that 

Sb{u) ^rhu- zu + 2al{2u'^ + 1), z > 0. (131) 

Taking into account the Jacobian of the transformation u = sinh s we arrive at 

ImZ^(i,m;a8) ~4v^^^e-(™-^)'/^^'^i+2'^i. (132) 

m — z 

Combining this with the fermionic integral we obtain for the spectral density 

[m — z) >y rh — 8al 

C. Edge Scaling 

In this section we study the Dirac spectrum near the gap aX z = Zg. We first consider 
z < Zg and then analyze z > Zg. 

1 . Edge Scaling for z < Zg 

We now expand Sbiu) near the edge of the spectrum. The second derivative vanishes and 
the first derivative is given by 

S',{Ug) =Zg-Z (134) 

resulting in the expansion 

Sb{u) = Sb{Ug) -{U- Ug){z - Zg) + ^(m " U gf S'^' {u g) (l35) 



with 



The saddle points are given by 



We will see that the minus sign corresponds to the leading saddle point of the imaginary 
part of the bosonic partition function whereas the negative sign corresponds to the leading 
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saddle point of the fermionic partition function as well as the real part of the bosonic 
partition function. Notice, that because of the supertrace, the fermionic and bosonic actions 
are each others inverses. The bosonic exponent at the saddle point for the imaginary part 
becomes 

S,{u) = S,{u,) - \^^^%r=^, (138) 
whereas the fermionic exponent at the saddle point reads 

Sdu) = -SJu„) - iMf^^ill!^ . (139) 

The second derivative at the saddle points is given by 

S"{U) = S"{Ug) + {U-Ug)S"\Ug), 

= iU-Ug)S"'iUg), (140) 

SO that 

Sbiu) = Sb{u) + ^{u-u)^{u-Ug)S"'{ug) 

= Sbiug) + i(n - - Zg)S"\ug)Y/\ (141) 

Therefore, integration over u gives an overall factor i for z < Zg for the saddle point u < Ug, 
whereas for the saddle point {u > Ug) the Gaussian integral is real. 
The spectral density is given by 

p,{z) = -Zf{z)^lmZ\z). (142) 
TT az 

Differentiating the pre-exponential factors gives subleading corrections. The leading order 
saddle point result for the spectral density is thus given by 

P5(5) = ^^e-i(^(^-^))^^Vv/=^. (143) 

— ^ [Ug) 

Note that the pre-exponential terms are not determined by a leading order mean field cal- 
culation. Comparing this result to the leading order asymptotic expansion of the universal 
result for the spectral density at the soft edge 

I g-4x3/V3 

pJx) = — (Ai'(x)^ - xAi^(x)) ~ , with x = (zq - z)/A (144) 

A TVX 

we find that the two results coincide if we make the identification 

A = {-S"'{ug)/2y/^ (145) 



42 



2. Edge Scaling for z > Zg 



For z > Zg the leading order saddle point result is determined by a pair of complex 
conjugate saddle points. Near the edge of the spectrum, they are given by Eq. f ll37p . 

In order to understand the mean field limit for z < Zg we first discuss the case of cig = 0. 
The saddle points are given by 

ZZ ITTL 

sinhf = ± — , coshf = =F — . (146) 

V -2^ — rh? yjz^ — rh? 

When the fermionic and bosonic saddle points have opposite sign, then the pre-exponential 
factor vanishes. Therefore this combination of the saddle points does not contribute to the 
partition function. Only one of the bosonic saddle points can be reached by deforming 
the integration contour and the relevant fermionic saddle point necessarily has the same 
sign. This is the way the supersymmetric method selects the fermionic saddle point and 
circumvents the failure of the fermionic replica trick. 

The pre-exponential factors add up to ^ z^ — which is canceled by the contributions 
from the Gaussian integrals about the saddle point (up to a factor of vr) resulting in a 
partition function that is correctly normalized. 



Next we consider the thermodynamic limit at fixed a. Then the 1 / yW—m^ singularity 



at the edge of the spectrum turns into a z — Zg singularity for which we expect an Airy 
like behavior. The saddle points for z close to Zg are given by 

(Kz - z)\^^'^ : 

sinhf-t = ± ? I — ^ I , coshf = y 1 + sinh f. (147) 



-S"\Ug) 

The structure of the saddle point approximation to the resolvent is given by 

G{z) = P{f+, f+)e^>'^'+^+^f^'+^ + P(f+, ^_)e^'>(^"+)+^/(-"-), (148) 

where the exponents are given in Eq. (!69|) and the prefactor is equal to 

(P(s,r) - ^) sinh r (149) 

with P{s, r) defined in Eq. (1701) . In the first term in Eq. f ll48p the exponents cancel, whereas 
in the second term they add up to 



Sbif+) + g.(f_) = I- (150) 
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When the saddle points are different the prefactor is suppressed by 1/m in agreement with 
the asymptotic expansion of the universal expression for the level density near the edge in 
terms of Airy functions. We have not worked out the asymptotic behavior for z > Zg but 
checked numerically that the exact expression for the spectral density is in agreement with 
the expression in terms of Airy function. The asymptotic domain is already reached when 
m ~ 40 and ~ 1- 
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Appendix C. THE SIGN OF AT SMALL MASS 



In this appendix we investigate in detail the sign of the partition function as a function 
of Ws, in the vicinity of m = 0. Throughout this appendix we will set z = 0. In particular 
we can show analytically that at m = the Nf = 2 flavor partition function is positive 
for all values of u only if Ws > 0. In contrast the partition function becomes negative for 
odd values of u if Wq < 0. For both signs of Ws we show that the Nf = 2 flavor partition 
function has a local minimum at m = 0: its first derivative with respect to rh vanishes, and 
its second derivative is positive. 

We begin by stating an equivalent integral form of the Nf = 1 flavor partition func- 
tion Eq. flMl) . by linearizing cos^(6') through a Gaussian integral and applying an integral 
representation for Bessel-/ functions: 

— e'^" exp \m cos(^) T 4a^ cos^(^)l 

e-" J, [m + 4a8x| ^ | J . (151) 

Here the sign +(— ) corresponds to positive (negative) Ws- Following Eq. (!93|) the A'^^ = 2 
flavor partition function is then given by 

At m = the integral in Eq. f ll5ip becomes doable. First, due to the parity of Bessel—/ 
functions, I^{—x) = {~yi^{x), the zero mass single flavor partition function vanishes for 
odd values of z/ (for both signs of Ws): 

Z';,=l\^\m = 0;±dl) = 0. (153) 

At even value of u we obtain for a single flavor 

■2k 

Z^^=i\{m = 0; ±dl) = { \ }lki2dl). (154) 

Thus for Nf = 2 at even u we have 

Z^=;j,{m = 0; ±dl) = { }Wdlf, (155) 

as only the first term in Eq. fll52p is non-vanishing. This is obviously real and positive being 
a complete square, for both signs of Ws- On the other hand when = 2A; -|- 1 the first term 
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vanishes in Eq. fll52p . and we have 



7u=2k+l 
'Nf=2 



{m = 0;±al) = {^ h{2al) h+,i2al) . (156) 



Only for the upper sign + the factor (— )2'«+i compensates the second minus sign. Knowing 
that Bessel-J is a positive function for positive arguments we get again a positive partition 
function for Ws > 0, but a negative one for Wg < 0. 

Next we investigate whether or not the point m = is a relative extremum. Using the 
following Bessel identity, 

lUx) = ^(/._i(a:) + /,+i(x)) (157) 

it is easy to see that 

-Z^,{m;±al) = Uz';{0;±al){Zr\0;±al) + Z'(^\0-,±al)) (158) 

', rh=0 L \ 



d_ 
dm 



Because of its parity this expression vanishes for both even and odd independent of the 
sign of W^. 

To see if this is a local minimum we compute the second derivative, given by 



^Z,^(m;±a^)|^^^ = ^ (2Zr(0; ±a^)^ - 2Zr^(0; ia^)^^ (0; ±a^) 

+Zri(0; ±alf - Zl~\Q- ±al)Z{+\Q- ±al) 
+Zri(0; ±alf - Zri(0; ±al)Z^,^\Q- ±al)). (159) 

For even v only the first line is non- vanishing due to parity, and we have 

ff. 1 , (_\2k 



Z]^='%m-±a,] 



m=0 



I }(4(2a^)2_4_^(2a2)4_,,(2a^)). (160) 



This is positive for both signs of 



h{xf -h-i{x)h+i{x)=2 f dtth{xtf > 0. (161) 

Jo 

Taking equation Eq. (11591) for odd z/, the first line vanishes due to parity and the second 
and third line do contribute to give: 

^Zr''^\m;±as)\^^^ = 1{ } (4(26^)^ - 4_i(2a^)4+i(2al)) 

+ 41 1 }(4+i(2«8)'-4(2aD4^.2(2a^)) . (162) 
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This is once more positive for both signs of Ws, because of the identity (116ip . In conclusion 
the second derivative is always positive and we have a relative minimum for small rh. 
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